Vertically coupled double quantum dots in magnetic fields 
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Ground-state and excited-state properties of vertically coupled double quantum dots are studied 
by exact diagonalization. Magic-number total angular momenta that minimize the total energy are 
found to reflect a crossover between electron configurations dominated by intra-layer correlation 
and ones dominated by inter-layer correlation. The position of the crossover is governed by the 
strength of the inter-layer electron tunneling and magnetic field. The magic numbers should have 
an observable effect on the far infra-red optical absorption spectrum, since Kohn's theorem does 
not hold when the confinement potential is different for two dots. This is indeed confirmed here 
from a numerical calculation that includes Landau level mixing. Our results take full account of 
the effect of spin degrees of freedom. A key feature is that the total spin, S, of the system and 
the magic-number angular momentum are intimately linked because of strong electron correlation. 
Thus S jumps hand in hand with the total angular momentum as the magnetic field is varied. One 
important consequence of this is that the spin blockade (an inhibition of single-electron tunneling) 
should occur in some magnetic field regions because of a spin selection rule. Owing to the flexibility 
arising from the presence of both intra-layer and inter-layer correlations, the spin blockade is easier 
to realize in double dots than in single dots. 
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I. INTRODUCTION 

Vertically coupled quantum dots have recently at- 
tracted much interest, since they open up the possibil- 
ity of manipulating electron configurations in a three- 
dimensional space. Quantum dots are generally thought 
of as 'artificial atoms', in which electrons are confined by 
an artificial electrostatic potential instead of being at- 
tracted to a nucleus. Numerical studies of few electron 
systems confined in single quantum dots have shown that 
the angular momentum and spin of the ground state in 
strong magnetic fields belong to a special series values 
called magic numbers. For single quantum dots the inter- 
val of magic number angular momenta has a one-to-one 
correspondence with the symmetry of charge correlation. 
This has been explained by an 'electron-molecule' picture 
by one of the present authors Q. Namely, the Coulomb 
repulsion forces the electrons confined in a dot to take a 
definite molecular configuration (triangle for three elec- 
trons, square for four, etc) when the electron correlation 
is strong [|l]J^]. If we use a similar picture for vertically 
coupled dots we can envisage the ground states in strong 
magnetic fields as three-dimensional electron molecules. 

On the experimental side, recent advances in nano- 
lithography and thin-film processing make it possible 
to fabricate vertically coupled multiple quantum dots, 
where two-dimensional (2D) electrons are confined within 
an area smaller than 1 /im across B. Observable ef- 



fects in the atom-like physics of dots have been detected 
from measurements of the tunneling current or capac- 
itance. In the low magnetic field regime (B < 2T), 
Tarucha et al have indeed found from capacitance spec- 
troscopy that Hund's rule and shell structures appear in 
the spectrum Q. As the magnetic field becomes stronger, 
however, the single-electron energy levels change to ones 
like those of a 2D harmonic oscillator to Landau levels, 
where levels with different angular momenta are degen- 
erate, and electron correlation should play an essential 
role. Ashoori et al and Wagner et al have looked at the 
addition spectrum (energies required for adding one ad- 
ditional electron) for single dots in magnetic fields, where 
a total-spin transition caused by the Coulomb interaction 
is expected 

This reminds us of the fractional quantum Hall ef- 
fect (FQHE) in the bulk 2D electron system, a mani- 
festation of strong electron correlation in high magnetic 
fields ||-|0|, where the total spin is a sensitive function 
of the density of electrons in the Landau levels. This is 
an effect of electron correlation, i.e., the way the elec- 
trons are correlated is strongly dominated by the total 
spin, while the Zeeman energy, which is few percent of 
the typical Coulomb energy, has only a minor effect. 

The fractional quantum Hall effect in double layers 
fil"[ has recently been under intense study. In this 
case the additional degree of freedom arising from the 
double layer (a pseudo-spin that labels the layers) en- 
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riches the physics. A central issue in these systems 
is the interplay of electron correlation and inter-layer 
electron tunneling. The competition between these ef- 
fects makes the quantum Hall state evolve continuously 
from a correlation-dominated (two-component) state to 
a tunneling-dominated (single-component) state within 
the quantum Hall regime. If we include spin degrees of 
freedom, FQH states specific to double-layer systems of 
electrons with certain specific spins indeed appear [ p^[ . 
Thus, vertically-coupled quantum dots are also intrigu- 
ing from the viewpoint of what happens if we laterally 
confine a double-layer FQH system. 

In the present paper we investigate the physics of dou- 
ble dots in magnetic fields. While we have studied the 
magic numbers and far infrared (FIR) absorption in our 
previous publications [j^ 18 1, we have assumed there that 
the system is fully spin polarized (except for our paper 
on the spin blockade, as recapitulated here in section |v|). 
We have also adopted the lowest Landau level approxi- 
mation there. These assumptions, however, are justified 
only for B — > oo limit. 

One of our most important findings is that the total 
spin of the ground state of the dots changes wildly as 
the magnetic field varies, where the ground state is spin- 
tmpolarized even for a magnetic field as large as B = 4 T 
in typical conditions. We show that this is a manifesta- 
tion of electron correlation rather than the spin Zeeman 
effect. Namely, the electron correlation energy is drasti- 
cally affected by the total spin state even when we ignore 
the Zeeman term, so that it is the former that determines 
the spin. A similar phenomenon occurs in single dots and 
is explained by the electron-molecule theory that takes 
care of spin quantum numbers Jl|||- The wild change in 
the spin can cause single-electron tunneling to be blocked 
in some magnetic field regimes due to spin selection rules 
as we have proposed recently [[l9| . Evidence that this ef- 
fect is caused by electron correlation is the occurrence of 
a re-entrant non-blocked region which is hard to explain 
by Zeeman-energy considerations. 

We have also investigated the effect of the difference 
in the curvature of the confinement potential in the two 
layers. This can be regarded as a pseudo-spin Zeeman en- 
ergy and a novel charge correlation caused by the pseudo- 
spin Zeeman energy is found to appear, which causes 
dependence of the magic numbers on the strength of 
the inter-dot electron tunneling and the layer separation. 
The upper branch in the absorption spectrum that ap- 
proaches the cyclotron frequency for large B has been 
obtained from a calculation that includes higher Lan- 
dau levels. We propose that the magic numbers should 
have an observable effect on both the higher- and lower- 
branches of the FIR spectrum. Specifically, while the 
generalized Kohn theorem states that the electron 

correlation cannot affect the optical absorption spectra in 
single dots with parabolic confinement, this is no longer 
the case with vertically coupled dots that have different 
confining potentials. We show that the absorption energy 
of the double dot should indeed exhibit discontinuities at 



the magnetic fields where the total angular momentum 
or total spin change from one magic number to another. 

Thus the purpose of this paper is two-fold: to include 
the spin degrees of freedom and to include the mixing 
of higher Landau levels, which enables us to calculate 
the real ground state. The magic numbers in double 
dots are found to depend on the strength of the inter- 
dot tunneling and this determines how the total spin (S) 
changes hand in hand with the total angular momentum 
(L), which can be distinct from the link between S and 
L in single dots. 

The organization of this paper is as follows: In Sec. || 
the theoretical model of a vertically coupled quantum dot 
is described. The physics of electron correlation in double 
dots and its relevance to the magic numbers is discussed 
in section Sec. Ill . The remaining sections deal with the 
experimental consequences of our results. All the results 
are summarized in Sec. VI 



II. VERTICALLY COUPLED QUANTUM DOTS 

Quantum dots are fabricated with a variety of tech- 
niques. The difference between these techniques lies 
in the way the lateral confinement potential of two- 
dimensional electron system is created. One method for 
obtaining lateral confinement is to etch away a semicon- 
ductor sample to obtain mesas which contain a laterally 
confined two-dimensional electron system. The dot struc- 
ture fabricated in this way is called a "deep-mesa-etched" 
quantum dot (Fig. [I] (a)). Another method is to deposit 
a metallic gate on top of a heterostructure. When biased 
negative the gate will deplete electrons from the region 
underneath it thereby creating a quantum dot (Fig. |l| 
(b)). Simple electrostatic considerations show that the 
potential felt by an electron in a mesa-etched sample is 
similar to that of a charged disk. The electrostatic po- 
tential of a charged disk can be found in terms of elliptic 
integrals and the bottom of the potential is well approx- 
imated by a parabola (2^j2^| . 

Recent advances in semiconductor fabrication tech- 
niques have enabled fabrication of double dots in vertical, 
triple-barrier structures on submicron scales. Fig. 0(a) 
actually depicts such a case: a gated field effect confined 
double dot and a deep-mesa-etched double dot with the 
source gate on top, which has been studied by Tarucha 
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As mentioned, we can assume that the confining po- 
tential of a quantum dot is parabolic. Then the single 
electron states are those of two dimensional harmonic 
oscillator. When a perpendicular magnetic field B is ap- 
plied, the Hamiltonian for a single electron in a quantum 
dot is given by 
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-7T H — m*uj n r H — q* unB<7 7 , 
2m* 2 2 y ^ B " 



(1) 



where m* is the effective mass, u>o the strength of the 
parabolic confinement potential, /i B is the Bohr magne- 
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ton, g* the effective g-factor, \a z the z-component of the 
spin of a single electron, and the canonical momentum 7r 
is given by 

7T = -ihV + eA, (2) 

where the vector potential A satisfies V x A = B. For 
parabolic confinement, the Hamiltonian is similar to that 
for a free electron, so that the eigenstates 4> n i, called 
Fock-Darwin states p^ , p5| , are given by 

/ n\ 
0n£(r) = y 27rA 2 2lWI(n+KI)! 

Here n is the radial quantum number, I the angular 
momentum quantum number, L„ a Laguerre polyno- 
mial, A = yjTi/m*Vl is the effective magnetic length, 
f2 = ywf + 4wq and lo c = eB/m* is the cyclotron fre- 
quency. The Fock-Darwin states are ring-shaped with 
a radius R ~ Ay2(2n + \i\ + 1). The eigenenergies are 
given by 

£ nto (0) - ~(2n + 1 + K|)M} - i«iw c + i<?*/i B Ba,. 

(4) 

For double quantum dot systems, the single-electron 
states are split into symmetric and antisymmetric states 
due to inter-layer tunneling (Fig. 2). The strength of 
the tunneling is characterized by the energy gap, A SAS , 
between symmetric and anti-symmetric states. Although 
A SAS depends (~ exponentially) on the layer separation 
(d), it also depends on the height of the barrier that 
separates the two layers. Experimentally, d and A SAS 
are independently adjustable. We can see that d controls 
the difference between the intra-layer interaction, oc 1/r, 
and the inter-layer interaction, oc l/(r 2 + d 2 ) 1 / 2 . 

The Hamiltonian for interacting electrons in a verti- 
cally coupled quantum dot, 

H = H + Ht+Hc, (5) 

comprises the single-electron part, 7io, and the Coulomb 
interaction, Hq- In second-quantized form with a Fock- 
Darwin basis we have 

= E EE £n/™4 ra c «'™> ( 6 ) 

n £ a a 
n t a 

(7) 
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(n 1 £ia 1 a 1 ,n 2 £2cr2a2\V{r 1 - r 2 )\n z £' i (T 3 ,a s ,ndiCF i a i ) 

XC n 1 £ 1 cr 1 a 1 C n 2 e2^2a2 Cn 'i- e 'i' J 'i a 'i Cn 3 e 3^30i 3 - (8) 

Here a is the real spin, a — ± is a pseudo-spin index 
specifying the two layers, , (c n e aa ) are creation (an- 
nihilation) operators and e is the dielectric constant of 
the host material. When the spin and pseudo-spin are 
included the energy of the single-particle states becomes 

Snlcra = Sn£cr(^a), (9) 

where E nl(y is given by eq.(f| and fl a = {uo 2 c + icJo a ) 1/2 
with woq the confining potential for the ath layer. 

Because the Coulomb interaction conserves the total 
angular momentum the many-body Hamiltonian is diag- 
onalizcd numerically in each sector of the total angular 
momentum space ^l| , ^6| -|29[| . 

III. MAGIC-NUMBER STATES 

One of the most dramatic features of interacting two 
dimensional electron systems confined in quantum dots 
subjected to magnetic fields appears in their energy spec- 
trum. Numerical studies of finite single-layer systems 
have shown that ground states only occur at certain to- 
tal angular momenta and total spins called "magic num- 
bers" . The reason why magic numbers appear is that the 
Coulomb part of the total energy is not a smooth func- 
tion of the total angular momentum, L, but has a series of 
cusps caused by electron correlation |H| . A larger L cor- 
responds to a larger spatial extent (~ y/L) of the charge 
density, which costs a higher confinement-potential en- 
ergy (Fock-Darwin energy) while the Coulomb repulsion 
tends to be reduced. What the magic numbers are all 
about is that the competition of these effects is not mono- 
tonic. Indeed, the magic numbers have a one-to-one cor- 
respondence with the symmetry of charge correlation, 
which has turned out to be the same as the symmetry 
of the classical equilibrium configuration of point charges 
for single dots. If we further take account of the electron 
spins, we have to consider, in addition to the Zeeman 
energy, the competition among the states having various 
total spins that determines how the electrons should be 
correlated. 

Magic number states also occur in double dots. For 
double dots, however, we have pseudo-spin as an addi- 
tional degree of freedom, which makes the sequence of 
magic numbers for double dots different from that of sin- 
gle dots. We first look at how the sequence of magic 
numbers is modified by the strength of inter-layer tun- 
neling. 
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A. Three-electron double dots 

For simplicity wc start with three fully-spin-polarized 
electrons in a double dot with the same confining poten- 
tial for both layers, while the effect of spin will be ex- 
amined later. We plot the ground-state energy (Fig. ||) 
calculated numerically as a function of the total angular 
momentum, L, at B = 15T for A S as = 0.2(0.6) meV 
in the upper (lower) panel. The confinement energy is 
Hluq = 3.0meV for both layers and the layer separation 
d = 20 nm. 

Magic numbers can be identified from the positions of 
downward cusps. For Asas = 0.2 meV we have a period 
of two (L — 5, 7, 9) for smaller L followed by a period of 
three (L — 9, 12, • ■ •). For a larger Asas = 0.6 meV, by 
contrast, the period is three for all L, as in the case of a 
single-dot containing three electrons. 

In the right panel of Fig. |3| we show the average charge 
density, p, which is axially symmetric. Since we assume 
the same confinement potential for the two dots, p is 
the same for the two dots. For L — 5 (solid line) the 
density against the lateral distance r from the center has 
a peak at the center along with a shoulder around 11 
nm (=1.7A, A: effective magnetic length is 6.54 nm). For 
L — 12 belonging to the period of three (long dashed 
line), on the other hand, the density is peaked at a finite 
r. 

In order to identify the mechanism for the change of pe- 
riod in the magic numbers, a better measure of the elec- 
tron correlation is the pair correlation function P(r, ro) 
which is proportional to the conditional probability of 
finding an electron at position r given that there is one 
at position ro. This is shown in Fig. |J. We fixed one elec- 
tron at the place where the charge density has a shoulder 
(or a maximum when the density is peaked at a finite r). 
We can immediately see that the ground-state electron 
configuration changes from the one dominated by inter- 
layer correlation to one dominated by intra-layer correla- 
tion. In particular, the correlation for L = 5 corresponds 
to a triangular "electron molecule" developed across the 
two layers, with one electron at the center of the lower 
layer while the other two are in the upper layer. In con- 
trast, the triangular form develops within each layer for 
L = 9,12. The L — 7 state is an intermediate one where 
the upper layer has a single peak but the lower layer has 
two peaks. 

An important difference between a single dot and a 
double dot is the way a magnetic field affects the plot of 
ground-state energy E against L. For single dots and 
strong fields such that states with n — and £ > 
are dominant, the Fock-Darwin energy in Eq. |] becomes 
h(p, — uj c )£/2 + Ml/2, The main effect of changing the 
magnetic field is then to change the coefficient of the 
term linear in I so that once we obtain an E — L plot for 
a specific value of B we can predict how the ground-state 
L will evolve with B. In other words, a series of cusps 
(magic numbers) in the E — L plot will sequentially be- 



come the true ground state as B is varied. For double 
dots a change in B affects not only the linear term but 
also the ratio of the effective magnetic length A and the 
layer separation d, or equivalently, the ratio of the intra- 
and inter-layer Coulomb interactions. Thus we have to 
check which of the magic number states can become true 
ground states. 

We have done this in Fig. ||. Here we have also in- 
cluded the effect of spin degrees of freedom. So what 
we have plotted in the figure is the total spin as well 
as the total angular momentum of the ground state of 
the three-electron double dot against magnetic field for 
A SAS = 0.2(0.6) meV in the left(right) panel. One can 
find that the magic- number states for [L = 5, S = 3/2) 
or (L = 7, S = 3/2) become the absolute ground state for 
B = 7.2 ~ 7.5T or B — 10.2 - 11.0 T, respectively for 
A SAS = 0.2 meV. On the other hand, the magic numbers 
L = 5, 7 are washed out when A SAS is increased to 0.6 
meV. 



B. Four-electron double dots 

The crossover between dominant intra-layer correla- 
tion and dominant inter-layer correlation also exists in 
four-electron systems. The total angular momentum and 
the total spin of the ground state of four electrons are 
shown in Fig. |[ The system is fully spin-polarized for B 
> 3.2T, where the magic-number orbital angular momen- 
tum evolves as L = 4, 6, 8, 10, 14. In this regime we can 
discuss the magic numbers without complications coming 
from real spins. 

The charge density in Fig. ^ is single-peaked for L up to 
14. However, unlike the three-electron case, this does not 
imply that the correlation is developed within each layer. 
If we look at the pair correlation function in Fig. |?], we 
can in fact see that the correlation is developed across 
the two layers. Since we have four electrons we have 
now a tetrahedral configuration (with two electrons in 
the upper layer, and the other two in the lower layer in 
a staggered position) for all the magic-number states. 

As the magnetic field increases, the period of magic 
numbers changes from two (L = 6 + 2 x integer) to four 
(L = 6 + 4 x integer) as shown in the left panel of Fig. ||. 
The absence of the magic number L = 12 is strong ev- 
idence for this change, because L = 12 belongs only to 
the sequence L = 6 + 2 x integer. The change of the 
period of the magic numbers corresponds to a crossover 
between dominant intra-layer correlation and dominant 
inter-layer correlation by analogy with three electron sys- 
tems, where the period of the magic numbers changes 
from two (L = 5,7) to three (L = 9,12) and L = 11 
is the missing magic number. The results for the pair 
correlation function in Fig. also show that inter-layer 
correlation develops as L increases. 
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C. Crossover between intra- and inter-layer 

correlations 

Intuitively the change in the correlation exemplified 
above for three and four electrons can be understood by 
considering the total energy. If the total angular mo- 
mentum L is decreased at fixed magnetic field the lateral 
spatial extent (oc \fT) of the wave function is compressed, 
so that the mean electron separation, a, becomes compa- 
rable with the vertical separation, d, of the layers. When 
the two layers are close enough, d <C a, the system is ef- 
fectively one-layer, so that all the electrons are correlated 
(with a square configuration for the four-electron case) 
within each layer. In the opposite case of separated lay- 
ers, d ^> a, the electrons in a, say, four-electron system 
should tend to minimize the energy by separating into 
two in the upper layer and another two in the lower end- 
ing up with a tetrahedral configuration across the layers. 
This is because the typical intra-layer Coulomb repulsion 
(oc 1/a) becomes significantly greater than that of the 
inter-layer Coulomb interaction (oc 1/ 'y 'a 2 + d 2 ), so that 
the electrons will avoid each other as far as possible by 
'sidetracking', i.e., by developing inter-layer correlations. 
Such staggered configurations have to involve the mixing 
of symmetric and antisymmetric states across the two 
dots via the inter-layer tunneling, and costs an energy 
Asas- This explains why the new magic numbers tend to 
emerge for smaller AgAs ■ Thus we can call the change in 
the nature of the electron correlation as a crossover from 
a dominant intra-layer correlation to a dominant inter- 
layer correlation. Further, if L is increased a becomes 
large and the inter- and intra-layer interactions become 
similar. Then the system is again effectively one layer 
and this explains the tendency for the new magic num- 
bers to occur for smaller L. These qualitative arguments 
apply for any value of the magnetic field but the mag- 
netic field determines the overall length scale and hence 
the L value of the true ground state. 

D. Effect of the pseudo-spin Zeeman energy 

When the two layers have different confinement poten- 
tials, the difference plays the role of a Zeeman energy for 
the pseudo-spin |i"7j]30[| , 

enia+ - e nta - = ~(2n + 1 + \l\)h(fl + - n_), (10) 

where f2± = (cu^+iuj 2 ^ 1 ^ 2 . The layer having the smaller 
confinement potential tends to accumulate charge and 
this causes a charge density imbalance between the two 
layers as shown in Fig. || for a ten percent difference in 
hu;o+ = 3.0 meV, TkJo- = 3.3 meV |j"7[ ]. For four-electron 
systems, the Zeeman effect for pseudo-spin leads to a 
novel magic number which corresponds to tetrahedral- 
like charge correlations |Q. 



Magic angular momenta and the total spin in this case 
are displayed in Fig. |[ We can see that the real spins are 
fully polarized for B > 4 T. Let us discuss this region in 
order to focus on the pseudo-spin. There the magic num- 
bers are L — 6, 9, 10, 12, 14, and we notice that the magic 
numbers L = 9, 12 do not appear in system without the 
pseudo-spin Zeeman energy. In particular, L = 9 be- 
longs to neither the sequence L — 6 + 2 x integer nor L — 
6 + 4 x integer. Instead it belongs to L = 6 + 3 x integer, 
which implies a three-fold symmetry in the charge corre- 
lation. 

The average charge densities in the ground states are 
shown in Fig. |^ for both upper and lower layers. The 
densities for L = 6, 9, 12, which belong to the sequence 
L = 6 + 3 x integer, are peaked at the center. If wc 
look at the pair correlation function in Fig. ||, the L = 9 
state is indeed seen to have a trigonal cone (three-by-one) 
configuration (while the L = 6 state may be thought of as 
a mixture of two- fold and trigonal configurations). The 
L = 12 state belongs to both L = 6 + 3 x integer and 
L = 6 + 2 x integer sequences, so that its correlation is a 
mixture of conical (three-by-one) and tetrahedral (two- 
by-two) symmetries. 

It would be interesting to compare our results with the 
phase diagram for the bulk bilayer FQH system obtained 
experimentally , or theoretically The picture 

in the bulk is that the quantum Hall state evolves, as 
AgAs is decreased, from a tunneling-dominated state to 
a correlation-dominated one. The former is a fully occu- 
pied symmetric state (one-component state), while the 
latter is an inter-layer correlated Laughlin's liquid, ^331 
(two-component) . 

To compare these with the dot states, however, 
is not straightforward. The bulk phase diagram is 
drawn against two dimensionless quantities, d/£s and 
AgAs/(e 2 /e^s), where l B = (h/m*^^ 1 / 2 is the mag- 
netic length. Because of the confining potential, the rel- 
evant length scale for dots becomes the effective mag- 
netic length A with A 2 = h/m*(uj 2 + 4W 2 ) 1 / 2 . With 
the parameters we have used A = (0.91 ~ 0.97)£b for 
B = (5 ~ 10)T. This yields e 2 /eA = 14.7 meV for 
B = 10T, so that A SA s/(e 2 /eA) = 0.01 - 0.04 for the 
double dots considered here. The Landau level filling, v, 
which is usually defined as 

v = N(N - 1)/2L 

for a dot having N electrons ranges from v = 3/5 for 
L = 5 to v = 1/4 for L = 12 for N = 3. Although there is 
some attempt at extending the bulk phase diagram from 
the usually studied v = 1 case to v = l/(odd integer) 
|l3| , wide variation in v and the fact that a change in B 
also changes A complicate the comparison. 

IV. FIR ABSORPTION 

Now we move on to the FIR absorption spectrum. In 
a single dot with a parabolic confinement potential, the 
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electron-electron interaction cannot affect the FIR ab- 
sorption. This is dictated by Kohn's theorem J33|, orig- 
inally conceived for translationally symmetric systems 
and later generalized to 2D dots with parabolic confine- 
ment in magnetic fields ■ Namely, long- wavelength 
electromagnetic radiations with electric field E couples to 
the dot via the dipole interaction, 

JV 

W' = ^eETi, (11) 

i=i 

which depends only on the center-of-mass coordinate. In 
a single dot with parabolic confinement the Hamiltonian 
separates into the center-of-mass and relative (interac- 
tion) parts even in the presence of magnetic fields, and 
the latter is irrelevant to optical transitions. Recent op- 
tical measurements of quantum dots indeed exhibit ab- 
sorption frequencies that are independent of the number 
of electrons and well fitted to the single-electron absorp- 
tion spectrum p4[ . 

In contrast, the separation of the center-of-mass and 
relative parts does not occur in vertically coupled dots 
having different confinement potentials even when they 
are both parabolic. This means that the Coulomb inter- 
action can affect FIR absorption spectra. 

To quantify the effect we have calculated the FIR ab- 
sorption spectrum of vertically coupled dots from the 
matrix element of the perturbation Hamiltonian, Ti.' , be- 
tween the ground state and all the excited states. Before 
discussing the results, we comment on the applicability of 
this approach to real systems. One important question is 
the nature of the electric field E in small samples. Several 
authors have questioned the relation between the applied 
electric field and the internal electric field in mesoscopic 
systems p5|-p7t with the general conclusion that depolar- 
ization effects are important. Therefore we would have to 
calculate the internal electric field to obtain the absolute 
value of the absorption coefficient. In addition, precise 
calculation of the absorption spectrum would require us 
to take account of other device properties that affect ab- 
sorption, such as finite thickness of the individual dots 
and deviations from a parabolic potential, about which 
scant information is available. We therefore make the 
reasonable assumption that the internal electric field is 
uniform and discuss only the absorption energy and the 
relative intensities of various transitions. This should be 
sufficient for our purpose of demonstrating that the FIR 
absorption of vertically coupled dots is affected by the 
electron correlation. 

The FIR absorption spectrum consists of two branches: 
the upper branch for inter-Landau- level transitions, 
which approaches to the cyclotron frequency in the B — > 
oo limit, and the lower one for intra-Landau level transi- 
tions. If we plot the position and intensity of the peaks in 
the FIR spectra along with the ground state total angu- 
lar momentum L and spin S (Figs. and |l|), we find a 
one-to-one correspondence between the magnetic fields 
at which the absorption line jumps and the magnetic 



fields at which the total angular momentum and/or the 
total spin of the ground state changes from one magic- 
number state to another. Thus the transitions should 
be directly observable in the FIR absorption spectrum. 
The figure also shows that the absorption intensity (oc 
square of the matrix element) is not monotonic. We have 
also displayed the FIR spectra for a four-electron sys- 
tem, which again exhibit similar jumps. The magnitudes 
of the jumps, obtained assuming a confining potential 
asymmetry of 10%, is a few tenths of a milli electron- 
volt. We believe they should be observable and the effect 
could probably be enhanced by making a double dot with 
greater potential asymmetry. 

V. SPIN BLOCKADE 

We now turn to the possibility of a spin-blockade. It 
has been known for some time that the Coulomb block- 
ade occurs in mesoscopic systems such as quantum dots. 
This is a combined effect of the discreteness of energy 
levels and the electron-electron interaction (charging en- 
ergy). Weimann et al then suggested that, if the total 
spin differs by more than 1/2 in the ground states for N 
and (N — 1) electrons, i.e., 

\S(N) - S(N - 1)| > i (12) 

the transport of an electron through the dot is blocked, 
and this should cause missing peaks in the conductance 
due to single electron tunneling at zero temperature. 
This is called the spin-blockade [ |38||39| ] and has been stud- 
ied theoretically in the case of no magnetic field, where 
we can apply the Hund's coupling picture Q. However, 
we cannot satisfy the spin-blockade condition for circu- 
lar quantum dots with parabolic confinement potential in 
this case, so that some modifications such as an anhar- 
monicity in the confinement potential |40j| or rectangular 
hard wall shape confinement potential ]38| have to be 
introduced. 

Now our idea here is to exploit the manifestation of 
the electron correlation in the total spin in the present 
system. We have in fact seen that the electron correla- 
tion dominates the total spin (S) of the ground state in 
a peculiar manner, as shown in figures ^, ^j, and ^, which 
comes from the fact that the magic values of L and S are 
linked. This is in a sense no surprise, since the total spin 
dictates, through Pauli's exclusion principle, the way in 
which the electrons can correlate, while the magic num- 
bers in L come from the electron correlation as one of the 
authors has shown from the electron molecule picture for 
single layer quantum dots. 

The correlation effect in fact happens when the typical 
Coulomb energy is much greater than the single-electron 
level spacing. Effects of electron correlation on the spin 
states via Pauli's principle are known to occur in ordinary 
correlated electron systems such as the Hubbard model, 
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but the present case is a peculiar manifestation in strong 
magnetic fields. 

Here we propose to utilize this electron correlation ef- 
fect to realize a spin blockade. For that we can look at our 
results for 5 to see if the spin-blockade condition, Eq. [l^, 
is satisfied in some regions of B. In addition, if there are 
low-lying excited states through which an electron can 
be transferred, the blockade will not occur except at the 
temperatures lower than the excitation energy. Thus we 
have to study the energies of excited states as well. We 
have plotted these along with |5(4) — 5(3) | for double 
dots that contain three or four electrons with the same 
parabolic confinement potential for both layer. The to- 
tal angular momentum and the total spin of the ground 
state are shown in the right panel of Fig. || for N = 3, 
and in the left panel of Fig. || for N — 4. 

The difference between the total spin of three and 
four electron systems, shown in the bottom panel of 
Fig. [l2|, indicates that the spin-blockade condition, 
1 5(4) - 5(3)| > |, is satisfied for 2.0 < B < 3.1 T and 
B = 3.7 T. In that region the excitation energies (top 
and middle panels of Fig. |l2|) for N — 3 and N = 4 
are of the order of 0.05meV, which is unfortunately only 
half the typical experimental resolution ~ O.lmeV. How- 
ever, it is possible to enhance the excitation energy by 
tuning hujQ, d, or A SAS . For instance, a double dot with 
TiLJa = 6.0 meV, d = 16 nm, and A SAS = 1.2 meV has 
the excitation energy of about 0.12 meV, which should 
be large enough to cause an experimentally observable 
spin-blockade between N = 2 and N = 3 |41f| . We also 
notice a level crossing between the second and the third 
excited states around B = 2.4T for N = 3, which should 
appear in the addition energy spectra. 

VI. CONCLUSION 

We have studied double quantum dot systems by using 
exact diagonalization including spin degrees of freedom 
and higher landau levels. We have found novel magic 
numbers specific to double dots and shown that they are 
related to the symmetry of charge correlations. In addi- 
tion, we have shown that we can change the charge corre- 
lation, total angular momentum, and total spin by vary- 
ing the strength of the magnetic field B. These changes 
have been shown to affect the optical absorption spec- 
trum and single electron tunneling. For the optical ab- 
sorption spectrum, we found that because of the break- 
down of Kohn's theorem, the spectra differ from that of 
non-interacting system with splitting and jumps in the 
optical absorption spectrum. The magnetic field where 
the spectrum has jumps is the field where the total angu- 
lar momentum and/or the total spin of the ground state 
shift from magic value to another. In the case of sin- 
gle electron tunneling, we have shown that the tunneling 
is blocked because of the spin selection rule — the spin 
blockade, in some magnetic field regions. 
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FIG. 1. A deep-mesa etched (gated) vertically coupled 
quantum dot is shown schematically in a(b). 

FIG. 2. Schematic single-electron wavefunctions in the 
direction normal to the layers in a double quantum dot. 

FIG. 3. Left: Ground-state energy against the total an- 
gular momentum, L, in vertically-coupled dots with three 
spin-polarized electrons for Asas = 0.2(0.6) meV at B = 15 
T in the upper (lower) panel. The confinement energy is 
Tiujo — 3.0 meV for both layers, and the layer separation is 
d — 20 nm. Arrows indicate the positions of the cusps. 
Right: Charge density (p) per electron for each layer for 
L — 5, 7 (belonging to period two) and L = 9, 12 (three) 
with Asas = 0.2 meV. 

FIG. 4. a) : Contour plot of intra-layer (upper panels) and 
inter-layer (lower panels) pair correlation functions, P(r, ro), 
for L = 5, 7, 9, and 12. One electron (filled circle) is fixed 
at the place where the charge density has a maximum in the 
upper layer. Thus the total charge is 1/2 (3/2) in the upper 
(lower), since we have 3/2 + 3/2 electrons while we fix one 
electron in the upper layer. An area with a linear dimension 
of 12A = 78.5nm (A: effective magnetic length) is displayed. 
The confinement energies, the layer separation, Asas, and 
B — 15 T are the same as in the previous figure, 
b): Bird's eye view plot of (a). A white sphere indicates the 
position of the fixed electron. 

FIG. 5. The total angular momentum L and the total spin 
S of a double dot containing three electrons with confinement 
energy ?kJo+ — fwo+ = 3.0 meV and layer separation d — 20 
nm. A SAS = 0.2(0.6) meV for the left(right) panel. 



FIG. 6. Left: The total angular momentum L and the 
total spin S of a double dot containing four electrons with 
confinement energy Tiuo+ = Tiujo- = 3.0 meV. The layer sep- 
aration is d — 20 nm and A SAS = 0.6 meV. 
Right: Charge density (p) per electron for each layer for typi- 
cal ground states in the fully spin-polarized region, B > 3.2T. 

FIG. 7. a) : Contour plot of intra-layer (upper panels) and 
inter-layer (lower panels) pair correlation functions, P(r, ro), 
are displayed for each of the magic- L states in the true ground 
state for the four-electron double dot shown in Fig. [] for ap- 
propriate values of B (in contrast to Fig. 4 where B is fixed). 
One electron (white sphere) is fixed in the upper layer at the 
place where the charge density has a maximum. Thus the to- 
tal charge is 1 (2) in the upper (lower) layer, since we have 2 
+ 2 electrons while we fix one electron in the upper layer. An 
area with the linear dimension of 12A (A: effective magnetic 
length) is displayed. The confinement energies, the layer sep- 
aration d and Asas, are the same as in Fig. ^| 
b): Bird's eye view plot of (a). A white sphere indicates the 
position of the fixed electron. 

FIG. 8. Left: The total angular momentum L and the 
total spin S of a double dot containing four electrons with 
confinement energy huio+ = 3.0 meV and hujo- = 3.3 meV. 
The layer separation is d = 20 nm and A SAS = 0.6 meV. 
Right: Average charge densities (p) per electron for upper 
(lower) layer in the ground states are displayed in the upper 
(lower) panel for L > 6, for which the spins are fully polarized. 

FIG. 9. Contour plot of intra-layer (upper panels) and 
inter-layer (lower) pair correlation functions, P(r, ro), for the 
same ground states as shown in Fig. ^[ One electron (white 
sphere) is fixed in the upper layer at the point where the 
charge density has a maximum. An area with the linear di- 
mension of 12A is displayed. The confinement energies, the 
layer separation d, Asas, and the magnetic field B are the 
same as in Fig. ^| 

b): Bird's eye view plot of (a). A white sphere indicates the 
position of the fixed electron. 

FIG. 10. FIR absorption spectrum (top panel), the total 
angular momentum (middle), and the total spin (bottom) of 
vertically coupled dots for N = 3 electrons. huio+ = 3.0 
meV, fkuo- = 3.3 meV, the layer separation d — 20 nm and 
Asas = 0.6 meV. The position of each filled circle gives the 
energy of the transition while the size of the circle represents 
the relative intensity of absorption. Vertical dashed lines are 
guides to the eye. 
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FIG. 11. FIR absorption spectrum (top panel), the total 
angular momentum (middle), and the total spin (bottom) of 
vertically coupled dots for N = 4 electrons is shown in the 
lcft(right) panel. Tkvo+ = 3.0 meV, fkvo- = 3.3 meV, the layer 
separation d — 20 nm and Asas = 0.6 meV. The position of 
each filled circle gives the energy of the transition while the 
size of the circle represents the relative intensity of absorption. 
Vertical dashed lines are guides to the eye. 

FIG. 12. Top(middle): Low-lying excitation energies for 
N = 4(iV = 3) double dots. 

Bottom: The absolute value of the difference, \S(4) — S(3)|, 
in the total spin for N = 4 and N = 3 double dots. 
hua + = Tiu)o- = 3.0 meV, the layer separation d = 20 nm and 
Asas = 0.6 meV. 
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